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(_^ ' Abstract. A Nash equihbrium has become important solution concept for analyzing the decision making in 

Cn I Game theory. In this paper, we consider the problem of computing Nash equilibria of a subclass of generic 

^.^i finite normal form games. We define rational payoff irrational equilibria games to be the games with all rational 

C^ ' payoffs and all irrational equilibria. We present a purely algebraic method for computing all Nash equilibria of 

these games that uses knowledge of Galois groups. Some results, showing properties of the class of games, and 

an example to show working of the method concludes the paper. 

o. 
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H ^ . 

rh \ 1 Introduction 

, ^, ' The problem of computing ah equihbria of class of integer payoff irrational equilibria(IPIE) games is stud- 
ied in [Ij. The approach in [Ij is as follows: They characterize equilibria of an IPIE game as solutions 

'T^ " to a system of polynomial equations called a game system QS. With a numerical method they compute a 

f^ — . ■ sample solution of the QS and follow it up by the Galois group action for computing all equilibria solutions. 

^D . Existing methods for computing equilibrium, such as the approach based on the homotopy continuation 
methods, given in ^2j,have the drawback of providing solutions via approximation. On the other hand the 
method in [3] is highly dependent on a chosen probability distribution. 

o: 

^D , With the motivation and game model presented in [I], in this work, we extend the approach of using Galois 
^7^ group for computing all equilibria of the IPIE games to a larger class of rational payoff irrational equilib- 

K*" ria(RPIE) games. The approach is a purely algebraic and provides a marginal decrease in computational 
^ , time compared to the method based only on Grobner basis [4] . The method computes exact equilibria and 

5^ \ does not depend on probability distributions. We further show that for a subclass of games it computes 
equilibria in closed formj^ 

The number of equilibria in a finite normal form game increases exponentially in the size of the game [2] . 
This means any method for computing all equilibria, one at a time, are bound to be exponential. In the light 
of this result structural relations between equilibria solutions of a game can help us compute all equilibria 
in a better way. In this paper, our aim is to present a pure algebraic method that computes all equilibria of 
RPIE games using knowledge of its sample equilibrium. Thus, answering partially an important question 
of establishing relations between all equilibria solutions, raised by McKelvey and McLennan [5]. 

To our knowledge, a method for computing equilibria of RPIE games with Galois groups has not been 
considered earlier. 



in 



cd 



^ cf. Proposition [31 



2 S CHATTERJI AND R GANDHI 

Organization of the article is as follows. Preliminaries and the game model is presented in Section [2j The 
algorithm for computing all equilibria of RPIE games and related results are presented in Section [3l Section 
m contains a worked out example of RPIE game. 

2 Preliminaries 

The subclass of games that we consider are generic and are known to have finitely many odd number of 
irrational equilibria [6] . Our convention is to write totally mixed real- irrational Nash equilibria as irrational 
Nash equilibria. 

Definition 1. A finite normal form game with 

— its all payoff values rational numbers 

— its all Nash equilibria irrational numbers 

is called an rational payoff irrational equilibria(RPIE) game. 

The characterization of equilibria as solutions to a system of polynomial equation and the underlying game 
model, that we consider in this paper, is based on [1]. 

Let T be an RPIE game with n = |A^| players. Each player i has /c, > 2 strategies, \Si\ = ki. We write 
/C^ = X]r=i ^i- ^jij2...?n denotes the payoff received by player i when each player adopts strategy jm for 
1 < jm < km and m = 1, . . . ,n. The probability that player i chooses strategy ji E {1, 2, . . . , ki} is denoted 
by X* . By Definition [1] an RPIE game has 



< 4 < 1. (1) 



Moreover, for each player i, 



Expected payoff for player i, 



K4 



E 4. = 1- (2) 



fcl k2 k„ 






a,; — > > ■■■ > Aj^j^^^^j^Xj^Xj^. . .Xj^ 



(3) 



In a Nash equilibrium, the following holds: 



fcl fe2 fci-1 ki^i k„ 

1 ™2 «-l i+1 

for every ji £ Si and for every i E {1, . . . , n}. 



"'^EE--- E E •••E^k- — ■■ .^xixi...x-x-^...x 



(4) 



With (JH), ^ and ([T]) we get a system of polynomial equation of the form. 

fcl fc2 f<^i-l ki + l kn 

ii=ii2=i ji_i=ij,+i=i j„=i 

for every ji E Si and i E { 1 , . . . , n} . (5) 
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For the purpose of computing all equilibria an RPIE game we shall consider the system of polynomial 
equations of the form ([5]) and call it a game system QS. We shall be specifically interested in the following 
field arithmetic for the QS representing RPIE games. 

Definition 2. Let K D F be a finite extension of a field F. Then the Galois group G = Gal{K/F) is the 
set 

G = {a : K ^ K\a is an automorphism, (T{a) = a for all a £ F}. 

We now consider the particular situation when K is an extension of F of the form F{a) where a is a root 
of a polynomial p{x) £ F[x]. It is known that the Galois group G of the extension K acts as a permutation 
group on the roots of the polynomial p{x). 

In this work our assumption is that these Galois groups are known, and we further utilize them for 
computing solutions of the polynomial system. It is important to note that for computing a Galois group 
of n! elements, computation of its 2n generators is sufficient [7J. It is a procedure not requiring exponential 
computation . And so the assumption, that the Galois groups are known, is feasible. The polynomials that 
we consider are multilinear - in indeterminate variables - over the field of rational numbers Q and generate 
irrational equilibria solutions. Thus the polynomial system induces non-trivial Galois groups. 

3 Equilibria of RPIE Games 

In this section, we outline the algorithm to compute all irrational equilibria of an RPIE game. We further 
show correctness of the algorithm and present other results. 

3.1 Method 

Before formally presenting our algorithm, we briefly discuss the approach and the underlying assumptions. 
We assume that we have an RPIE game T. As described in Section [2l we can derive a system of polynomial 
equations QS whose solutions include all the Nash equilibria of the game T. Since some of the solutions need 
not be Nash equilibria, our algorithm rejects these unwanted solutions using various different mechanisms. 
Due to Bernstein's theorem [8], we get an upper bound on the number of solutions a polynomial system 
can have. Bounds on the number of equilibria solutions - out of all polynomial solutions ~ are available in 
[9|10j . An upper bound on the number of solutions and bounds on the number of equilibria solutions give a 
bound on the number of solutions to be computed and the number of non-equilibria solutions to be rejected. 

In the initial phase of our method, Buchberger's algorithm is called to derive a univariate polynomial in the 
Grobner basis(GB) of the ^5o Since the game is RPIE, Nash's theorem [12] guarantees that the univariate 
polynomial has at least one irrational root. A root of the univariate polynomial is computed and substituted 
in the triangular form of a Grobner basis to find a univariate polynomial in other indeterminate variables. 
We repeat this procedure at most /C^ — n times, and at the end of it we have an irrational solution of the 
QS, a sample solution. 

We denote a Galois group of the irreducible part of a univariate polynomial in GB of QS by G. We assume 
G's are known. 



■^ Recall that the system of polynomial equations QS over complex number field has finitely many solutions. This finite 
variety of the QS (or equivalently zero-dimensional ideal X of the QS) guarantees a univariate polynomial in its Grobner 



basis [TT] , 
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In the next phase, we apply the transitive Galois group action corresponding to each indeterminate variable 
and find Galois-orbits to determine all irrational solutions of the QS. 

The final phase consist of determining all non-equilibria solutions and rejecting them. For this we use the 
Nash equilibrium verification algorithm in [13]. Following is an outline of the algorithm to compute all 
equilibria of an RPIE game with Galois groups. 



Algorithm 3.1 Computing All Nash Equilibria of an RPIE game. 
Input: An RPIE game, Galois groups. 
Output: All equilibria of the input RPIE game. 

1: /3 = (^1, /32, . . . , /3/C+)- {Initialize an empty tuple to store a sample solution of the QS}. 

2: Characterize all the Nash Equilibria of the input game as solutions to the GS. 

3: Call Algorithm 13.21 with QS for computing a sample equilibrium of the input RPIE game. 

4: Call the Galois group action Algorithm 13.31 with the sample solution tuple saved in /3. 

5: Save output of the Algorithm 13.31 in X. 

6: Reject non-equilibria solutions of the GS from X using verification algorithm in [13] or criteria ([T]) and 



Algorithm 3.2 Computation of a sample solution. 

Input: QS of the input game. 

Output: A sample solution /3 = (/3i, /32, . . . , /3x;+) of the input game. 



1: With Buchberger's Algorithm on QS, compute triangular form of GB. 

2: while one sample solution /3 of the QS is not constructed do 

3: Compute a root a of univariate polynomial - of some indeterminate variable Xj - generated in Step 

3. 
4: if a G Q then 
5: Reject a and go to Step 3. 

6: else 

7: Save a in /3 at location /3j. 

8: end if 
9: Substitute the root f3i in /3 into QS and compute a new triangular form with one less indeterminate 

variable. 
10: end w^hile 



The group action algorithm is given in [T]. We reproduce it here for the discussion to be complete. Algorithm 
13.31 computes group action by transitive Galois groups. The action is computed for each indeterminate 
variable Xi by considering it over each coordinate root in the tuple /?. The action generates Galois-conjugates 
of the roots that are further saved in as solution tuples in the set X. 

Traditional approach, given in [3], for computing solutions of system of polynomial equations using GB 
calls the Buchberger's algorithm for computing a triangular form. The triangular form provides a univariate 
polynomial in one indeterminate variable. Each root of the univariate polynomial is then substituted back 
in the triangular form to compute corresponding solution tuple, requiring multiple substitutions and fac- 
torizations. Algorithm 13.21 on the other hand, invokes Buchberger's algorithm exactly once. The Algorithm 
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Algorithm 3.3 Computing orbit of a Galois Group Action. 

Input: A sample solution /3 of the GS. 

Output: All the conjugate solutions of the input sample solution in set X. 



6 

7 

8: 

9 
10 
11 



Initialize the processed-elements list X and unprocessed-elements list U as X = U = {/3}. 
while U is not empty do 

Let u = {ui,U2, . . . , Uj(^+) be the first element of U. Delete u from U. 

for each i and j, g^j in Galois group Gi and Ui £ u. do 

9) 



Compute the transitive Galois group action n! 

1 2 K"^ 

if P' ^X then 

X = XLI {13'} andU = UU {jS'}. 
end if 



end for 
end w^hile 



computes a sample solution tuple corresponding to first irrational root of the univariate polynomial. 
Rest of the solutions are then generated by polynomial time group action, requiring no further substitutions 
and factorizations. 

It is important to note that, due to [12] and the fact that the input game has all irrational equilibria, we 
are guaranteed to get one solution of the GS in /? and so the Algorithm 13.11 reaches Step 4, every time. 
Next it calls the Algorithm 13.31 for computing polynomial time Galois group action over available sample 
solution in the /3. In the Algorithm 13.31 all other conjugate roots are computed with their known Galois 
groups Gi. 

Moreover, finite group action on finite variety guarantees that the Algorithm 13.31 reaches Step 13. At the 
end of Step 13, Algorithm 13.31 generates solutions of polynomial system QS in X, all of which may not be 
Nash equilibria. We use polynomial time algorithm, suggested in [13], to reject the non-equilibria solutions. 



For checking the condition in Step 4 of the Algorithm 13.21 i.e., for deciding q € Q we consider the following 
approach. Compute an approximation of a with suitable numerical method. With the KLL Algorithm 
|14j compute the minimal polynomial for a. Check irreducibility of the minimal polynomial over Q with 
algorithm in |15j . The irreducibility test decides whether q G Q or not. 



The Algorithm 13.21 computes a sample solution of the GS. Various approaches for computing a sample 
equilibrium of a game are discussed in |5]. The approach that we consider of using Grobner basis, though 
not new, is developed independently. Contrary to the other approaches based on Grobner basis, our ap- 
proach, given in Algorithm 13.21 differs primarily in Step 4. 

The Algorithm 13.11 shows a method for computing solutions of a system of polynomial equations without 
having to factorize the system every time. 

3.2 Comparison with the Algorithm in [1] 

The Algorithm in [T] calls MVNRM for computing an approximate solution of the GS. MVNRM is followed 
by the KLL algorithm for converting an approximate root in to algebraic form. The solution is then tested 
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for irrationality. In case otherwise, a new approximate solution is computed. 

On the other hand, the Algorithm 13.11 uses Buchberger's algorithm for computing a Grobner basis(GB). 
This converts the GS in triangular form, requiring doubly exponential time in /C"*". A univariate polynomial 
in the triangular form is factorized to find an irrational root. The root is substituted back in the GB for 
finding a complete irrational solution. This requires rejection of rational roots, if any. 

The difference in both the methods is: the Grobner basis method is purely algebraic and allows lot more 
structure to manipulate QS. While the numerical MVNRM provides lot more efficiency, disclosing a little 
about the structure of the polynomials, and interrelations of its solutions. 

3.3 Results 

The Algorithm 13. II computes all equilibria of RPIE games with n > 3 players. Following results show, why 
it can not work for games with n = 2 players. 

Proposition 1. A bimatrix game with all rational payoff values has all rational equilibria. 

Proof. For computing all equilibria of a bimatrix game, system of linear equations are sufficient [16] . That 
is, if all the game payoff values are defined over field F, then all of its equilibria - characterized as solutions 
to system of polynomial equations - can be found in the field F only. 

D 
Following is an immediate corollary to the result above. 

Corollary 1. Algorithm \3. 1\ can not be used to compute equilibria of a bimatrix game defined over arbitrary 
field. 

Proof. From Proposition [H it is clear that the class of two player games do not produce field extensions. 
This means that corresponding Galois groups are trivial, i.e., Gal{¥/F) = {e}. This proves the claim. D 

Corollary 2. The class of RPIE games is empty for n = 2 players. 

Proof. Follows from Definitions [H [2] and Proposition [TJ D 

Next result shows correctness of the Algorithm 13.11 

Proposition 2. The Algorithm \3.1\ to compute equilibria of the class of RPIE games works, i.e., the 
Alaorithm \3.1\ generates all irrational equilibria and no other solutions at termination. 

Proof. An input RPIE game T with n > 3 players is characterized as the QS of form ([5]). Polynomial 
system comes from the inequality on expected payoffs and payoffs at pure strategies. This causes the QS 
to have more solutions then just the equilibria. 

In the first phase. Algorithm 13.11 calls the Algorithm 13.21 The algorithm 13.21 computes a sample solution /3 
by first building a Grobner basis for the QS using Buchberger's algorithm. Buchberger's algorithm termi- 
nates with triangular form analogous to echelon form in the linear case. 



Since the game is known to be RPIE and rational solutions of the QS are rejected by the Algorithm 13.21 
the sample solution /3 must have all irrational coordinates. Consequently, each coordinate /3j of the sample 
solution 13 results in an algebraic extenion iT = Q(/3j) of Q with finite Galois group Gj = Ga/(]K/Q). Since 
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the group action of Gi is transitive, it generates all irrational solutions of the GS. 

Zero-dimensional ideal of GS guarantees that the group action terminates. This enables the Algorithm 13. II 
to reach Step 5, every time there is a RPIE game T as its input. The algorithm generates solutions of the 
GS that contain all the equilibria of the game T. 

Finally, Algorithm 13.11 rejects solutions of the GS which are not Nash equilibria. Since the set of Nash 
equilibria is known to be non-empty, set X contains all and only the Nash equilibria solutions of the RPIE 
game T. D 

With the available finite precision technology for representing a number in computer memory, the problem 
of storing an irrational equilibria is important [TT]. We show that the issue for RPIE games can be resolved 
as follows. 

Proposition 3. If univariate polynomials in ideal I of GS of an RPIE game has solvable Galois group. 
Then the Alaorithm \3.1\ computes Nash equilibria of the game in closed form. 

Proof. Due to Galois correspondence, we know, if a polynomial has solvable Galois group, all its roots can 
be computed by radicals. If each univariate polynomial in ideal I of GS of an RPIE game has solvable 
Galois groupp then roots of this set of polynomials can be computed using radicals. This gives all solutions 
in closed form, a subset is Nash equilibria of the game. D 

It is known that all abelian groups, groups of order < 60, group with odd order (Feit-Thompson theorem) 
and groups of order p"'q are solvable, where p and q are prime. Moreover, some non-abelian group may 
also be solvable [19]. This suggests that Proposition [3] is applicable for a substantial number of games. The 
equilibria of games with non-solvable Galois groups can be obtained in algebraic form by first computing 
equilibria solutions numerically, and then constructing minimal polynomial for each of the numerical values 
with the algorithm in |14j . 



4 Example 

In this section, we show working of Algorithm 13.11 with an example of 3 players 2 strategy RPIE game. 
With the Membership algorithm in [1], we verify that the game, given in Table [1] is RPIE. 





A 


B 




A 


B 


a 


6,-1, 4 


0, 9, 


a 
b 


2, 0, 


0, 9/2, 


b 


0, 3/2, 


2, 0, 


0, 27/2, 


4, 0, 6 



Table 1. Payoff table of a 3 player 2 strategy RPIE game. Player 1 and 2's strategies are indicated by a, b 
and A, B respectively. Player 3's strategies are 1 and 2. Entry in each cell of payoff table indicates player 
1, 2 and 3's respective payoff for their respective strategies. 



We let X = x\,y = x^, z = x^he the first strategy of players 1, 2 and 3 respectively. Probability that players 
will choose their second strategy is 1 — x,l — y and 1 — z respectively. After characterizing Nash equilibria 
as the solutions of a GS of form ([5]), we apply Buchberger's algorithm, with x ~< y ~< z, to compute a 



Using polynomial time Landau- Miller Test 118 1 it can be verified whether given polynomial has solvable Galois group or not. 
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univariate polynomial 5x^ — 27a; — 27 = in a GB. The univariate polynomial has x = ^(—9 it vl41) as 
its two roots, and is irreducible over Q. Its Galois group ({id,conjugate}) is isomorphic to Z2. 

Substituting x = ^(—9 — \/T41) in the triangular form of the GB and solving for univariate polynomials 
in y and z we get, y = 2^(3 — y/l41);z = ]j(14 + \/T41), a sample solution. Galois group of irreducible 
polynomials of polynomial system are known a prior (isomorphic to Z2 for each variable x,y,z over Q). 
Rest of the polynomial solutions can be obtained by computing Galois-orbits of the sample solution. 

In this example it is simple to observe the criterion ([T]) and ([2]) for deciding what solution tuple will form 
the equilibriaci Accepting values between and 1 we get x = ^(—9 + \/T4l);y = ^(3 + \/T41);z = 
jj(14 — -v/T41). This solution constitutes unique Nash equilibrium of the RPIE game depicted in Tabled) 
Note that, Galois group for the game system for this game is solvable, and so, all the equilibria computed 
are in closed form. 



5 Computational Complexity 

A Grobner basis can be computed in doubly exponential time in the size of /C"*". A Grobner basis contains 
polynomials in triangular form, and we are interested in the equilibria points with irrational values. An 
advantage of the triangular form is that at every stage of the substitution, unwanted solutions can be 
filtered out. For further details of computational complexity of finding Grobner basis see [20j . 



The group action Algorithm 13. 3^ in the worst case, requires action of each of the Galois group generator 
g' £ G' 'Z G to each element of the set of roots jCLj. This gives worst case time 0(|G"| • \X\). If a univariate 
polynomial has n roots then \G'\ is linear in n |7|21j . while |^| is polynomial in n. 

The algorithm for computing Nash equilibria with Grobner basis, given in [4j, substitutes all the roots in 
triangular form and solves univariate polynomial for each substitution. If Galois group is known for the 
polynomials then our approach computes solutions with relatively simple and efficient group action. The 
feature in our algorithm exploits information available in a sample solution. Our approach performs better 
then algorithm based only on Grobner basis. If each univariate in Grobner bases of the game system has 
di distinct roots {i € {!,... ,/C"'"}), then method for computing Nash equilibria with only Grobner basis 
takes Ilidi substitutions and factorization. On the other hand, in our approach, after computing a sample 
solution, no further substitution or factorization is required. 



6 Conclusion 

We defined a class of rational payoff irrational equilibria games and proposed a purely algebraic method for 
computing its all equilibria. The method that we presented throws more light on structure of equilibria and 
shows a way to effectively utilizes knowledge of a sample equilibrium. We discussed results that show why 
our algorithm can not be used for any bimatrix games, and how it can be utilized to compute equilibria 
in closed form. Working of our algorithm was demonstrated with an example of 3 players 2 strategy RPIE 
game. For the class of games that we consider, the algorithm we suggest is new. It would be interesting to 
see working of our algorithm on graphical forms and extensive form games, with necessary modifications. 



For larger system the Nash equihbrium verification algorithm comes handy. 
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